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We study electrical and thermoelectric transport properties of the surface state of the topological 
insulator and graphene in the presence of randomly distributed impurities. For finite impurity 
strength, the dependence of the transport coefficients as a function of gate voltage, magnetic field and 
impurity potential, are obtained numerically. In the limit of zero impurities (clean limit), analytic 
results for the peak values of the magneto-oscillations in thermopower are derived. Analogous with 
the conventional two dimensional electron gas, the peak values are universal in the clean limit. 
Unlike graphene, in topological insulators the coupling of the electron spin to its momentum leads 
to a dependence of the transport coefficients on the gyromagnetic ratio {g). We compare our results 
with data on graphene and identify unique signatures expected in topological insulators due to the 
magnetoelectric coupling. 

PACS numbers: 75.47.-m, 72.80.Vp 
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I. INTRODUCTION 

Helical Dirac fermions, massless relativistic charged 
particles with spin locked to their linear momentum, are 
proposed to exist on the surface of three dimensional 
topological insulators (TIs)i"— and later confirmed in 
experiments^!^ on samples like Bi2Se3 and Bio.gSbo.i . For 
compounds such as HgTe and Bi2Se3, surfaces with a sin- 
gle Dirac cone have been found^i^. While the low energy 
spectrum is very similar to that of graphene, the differ- 
ence in the microscopic origin of the phenomena points to 
subtle differences. For example, the presence of a single 
Dirac cone in TIs, as opposed to two doubly degener- 
ate ones in graphene, leads to significant reduction in 
phase space at low energies. Furthermore, the coupling 
of physical spin, rather than pseudospin (related to sub- 
lattice symmetry), in TIs suggests that the response in 
magnetic fields will be markedly different. 

The key property of the surface states of TIs is the 
magnetoelectric effect. Unfortunately, the presence of 
bulk carriers has made it hard to observe it in trans- 
port measurements^rJ^. One way to overcome this is to 
identify unique signatures, in transport, of the surface 
states. In particular the anomalous dependence of trans- 
port coefficients on the gyromagnetic ratio (g-factor) can 
be exploited to this end. The advantage of this approach 
is that the anomalous contribution can be enhanced by 
applying an in plane electric field^"^. In this article we 
study the thermopower and magnetotransport proper- 
ties of these surface Dirac fermions in the presence of 
randomly distributed impurities. 

In the conventional two dimensional electron gas under 
a quantizing magnetic field, the diagonal thermopower, 
Sxx, of the clean sample at low temperatures shows a 
series of peaks near the Landau levels N . The peak val- 
ues are \D.2kB / e{N + 1/2) independent of magnetic field 
strengtbi^. One of our goals is to show how these uni- 
versal features are modified in graphene and TIs. 

For electrons exhibiting a Dirac spectrum, such as 



graphene, this height of Sxx is expected to vary as 
ln2kB /eN , as required by the Berry phase effects lead- 
ing to a i integer shift in the Landau level indesi^ii^. 
We verify this result for |iV| > 1 and discuss the sin- 
gular behavior near iV = 0. We use self consistent 
Born approximation (SCBA) to consider the effect of 
short ranged, randomly distributed impurities on trans- 
port and magnetotransport^^'^* of the Dirac fermions 
and compute diagonal and off diagonal thermopower nu- 
merically. Similar to the results in conventional 2D met- 
als, the height of Sxx at low temperature does not have 
universal value after the inclusion of impurities^'*. 

We also show that for intermediate impurity strength, 
not considered in previous studies which have focused 
on the unitarity limitii^— , the particle hole symmetry 
of the original Dirac spectrum is broken under SCBA 
while for weak and strong impurity strength the parti- 
cle hole symmetry is restored. Another source of scat- 
tering in graphene is charged Coulomb impurities and 
its effect on conductivity and thermopower, both within 
perturbation-^^'^* and SCBA*^'^^, have been extensively 



studied. Comparison to available datap 
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suggests 



that a single scattering mechanism cannot account for all 
observed features^ -. Our formulation for TIs, with the g- 
factor set to zero, is similar to the case of graphene. Thus, 
in addition to exploring the difference in the transport 
properties induced by the spin orbit coupling, we can also 
compare our prediction with the available data. Quali- 
tative agreement with thermopower data is obtained for 
all but the zeroth Landau level. The failure is related to 
the underestimation of the longitudinal conductivity. 

The important new feature in topological insulators, 
as compared to graphene, is the spin-orbit coupling. In 
addition to orbital quantization, an external magnetic 
field couples to the momentum via its Zeeman coupling 
with the spin. Such an interaction reveals itself in novel 
signatures in transport unique to topological insulators. 
For example, the universal amplitudes in thermopower 
depend on the gyromagnetic ratio (g), but no splitting of 
peaks occurs as the surface states are derived from a sin- 
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gle Dirac cone. While the dimensional coupling constant 
a = g^fJ'^/vpe is quite small for HgTe (of order 10~'^/T) 
to yield a measurable signature, we derive the expected 
behavior in a model system where a is enhanced by an 
order of magnitude. The enhancement is achievable in 
principle by applying an inplane electric field^'^. 

The paper is organized as follows: in SeclUwe describe 
the general Hamiltonian and T-matrix formulation used 
to perturbatively compute the effects of impurities. In 
Sec. Uni we present our results, within linear response, 
for conductivity and the diagonal thermopower S^x in 
the absence of magnetic fields. In Sec. IIVI we consider 
the linear response conductivity and the thermopower 
Sxx and Sxy in the presence of the magnetic field. We 
compare our numerical results with those obtained in the 
graphene experiments and other theoretical results. Con- 
clusion are summarized in Sec. fVl 



where a;„ = 27r(n + 1/2)//? is the fermionic Matsubara 
frequency and 

The T matrix T(a;„) is 

tk) = ^"4i-^I]gO(^,c<.„)]"' (4) 

k 

In the following two sections we discuss the nature of 
transport without and with the magnetic field. We use 
the above T-Matrix formulation to compute the Greens 
function. Transport coefficients are obtained within lin- 
ear response using the appropriate current-current corre- 
lation function. 



II. HAMILTONIAN AND T MATRIX 
FORMULATION 

The surface state of the electrons on a topological in- 
sulator under electric and magnetic fields is— 

iJo = J drxl:\r)[vFd ■ TT — III — giiBS ■ B — eEx\'4){r){l) 

where Wi? is the Fermi velocity {h is restored in computing 
current current correlation), a = axX + Uyij denotes the 
Pauli matrices in spin space, / is the 2x2 identity matrix, 
i!{r) = ('0-1-, -04,)^ is the annihilation operator for Dirac 
spinor (T denotes the transpose of a row vector) and tt — 
— iV 4- eA is the canonical momenta. Electric field is 
assumed to be pointed in x-direction and magnetic field 
is perpendicular to the plane (z-direction). The coupling 
to local impurities of strength V is modeled as 



Hi, 



dr'^ V5{r- fi)tl)l{r)'il)s{r) 



(2) 



where s denotes spin degree of freedom of electron opera- 
tor ^{r) and sum over indicates summing all impurity 
positions. The full Hamiltonian of the system is given by 
H = Hq + Himp- In the dilute impurity limit we compute 
transport coefficients to linear order in impurity density 
rii = Ni/N. The finite temperature Green's functions of 
the electrons are 



where k — {kx,ky) is the two dimensional momentum of 
the surface band, s and s' are the spin indices, and 

G,,v(fc,P,r) = -(T0,,^(r)Vj;,.'(O)) 

with T being the time ordering operator. To linear order 
in Ui, the Greens function satisfy 

W„) = G"(fc, iUn) + G°(fc, W„)T(w„)G(fc, iOn) (3) 



III. ZERO MAGNETIC FIELD 

In zero magnetic field and i? ~ the Hamiltonian for 
the helical metal without impurities reduces to 



Hq = / drip'lrjvF 



fi kx 



The eigenfunction of the Schordinger equation 



Ho J drFj:^^{r)cl^\0) = {jvpk ~ fi) J drFf: Jr)ci _^\0) 
in the energy eigenstate is given by J dfF^ ^{f)ct |0). 



Here |0) is the vacuum, cl denotes spinless fermion op- 

erator, |fc| — fc , and the spin part is described by the two 
component spinor function ^{r) is 



1 



/2A 



7^ 



(5) 



with 7 = +/ — 1 denoting conduction/valence bands, 
(f>{k) = tan~^{ky/kx), and A is the area of the system. 
Following the full self consistent Born approximationii 
we have the self energy of electron given by 



1 - Iv T.k^'yl-'^ + M - jvpk - S(w)] 



The derivation of Eq.® is in the Appendix|Xl This self 
energy is solved numerically as shown in FiglT] The main 
feature of this self energy is that for moderate impurity 
strength {V — lOeT^) the particle hole symmetry is not 
preserved, while for strong {V — lO'^ey) or weak {V — 
10~^eV) impurity interaction the particle hole symmetry 
is restored. The scale is set by the bandwidth which in 
our calculations is 6eV. 

The physical quantity we are interested in is the trans- 
port properties of surface states. In the Kubo formula- 
tion the linear response conductivity is related to current 
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FIG. 1: Real part (left) and imaginary part (right) of the 
self energy plot. D — 3eV, riimp = 10~^, and V = WeV for 
blue lines and V = lO^eV for dashed purple lines. The inset 
of right figure shows imaginary part of self energy for weak 
impurity potential V = 10"^ eV. 



current correlation. Written in Matsubara formulation 
the particle current current correlator ia^ 



contributions to the conductivity comes from non cross- 
ing diagrams. The momentum integral k is cut off by 
linear spectrum energy boundary D with kcVp = D. 
This upper cutoff is also relatedii to unit cell size V by 

The zero frequency thermal response function L^j is 



related to L]] hy^ 



= lim 



dee- 



and the thermopower Sij is 



1 

raj 



(12) 



. L/(c^,r = 0,e + /i), 



(13) 



-iT 



h{iu})NV Jo 



h/3 



dTe''^^MT)Jj{0)) (7) 



Here NV = A is the area of the system, with N being 
the number of unit cells and V denoting unit cell area. 
Ljj is related to the linear response charge conductivity 



an via 



(8) 



The particle current operator in k direction in the 
Heisenberg representation is 



Jk = -77- = i[H, rk\ = -VpfTk 
at 



(9) 



Note that the electron current is related to the spin of 
the electron and hence sensitive to the Zeeman coupling. 
The corresponding particle current operator in the energy 
eigenbasis, Eq.([S]), is 



k,f K.7 



Jy = i 



.VF 



^(^e^*W_^e-^(^^))4^^c,,- (10) 



Within the limit of dilute random impurities we use the 
SCBA which is valid as long as kpl ^ 1 (with I as elec- 
tron's mean free path)^''. The Kubo formula for the real 
part of particle current current correlation is then given 
as 



— [nF(e) - npit + uj)] 



X 2^ 3G(fc, e iO+)3G(fc, e w + iO+) (11) 



with np{e) — 1/ {e^*^ + 1) (Fermi Dirac distribution func- 
tion) and Ll.y = 0. Quantum effects such as weak local- 
izations are not considered in this article as the major 



In zero magnetic field only diagonal component of Sij is 
nonzero. For a given self energy I](e) = 3?S(e) + i3I](e) 
the DC conductivity is 



/ , dnp{e — fj.) 



(14) 



3E(e) 



X / d{vFk){vFk)(- 

»S(e) 



— f , dnpie — /i) 



2'Kh 



de- 



de 



'^(Cj IT'imp) 



The explicit form of K(e, V, riimp) is given in the Appendix 
IbI Notice that for ~ the low temperature conduc- 
tivity is proportional to ~ 0, V, riimp) — 1- Thus the 
low temperature conductivity at Dirac node has an uni- 
versal value e^/27r/i. Close to zero temperature Eq.([8]) 
to Eq. p^ lead to the generalized Mott formula for the 
thermopower— 



3e 



(15) 



For a clean surface state we take V and nimp — and 
thus away from the half filling we have D ^ /i— 5RS(/i) ^ 
5S(/i) ~ 0. Near zero temperature in the clean limit the 
thermopower takes the form 



Sx 



kB ( 1 - a^3fiS(Ai) 



3e 



{^l-m:{^x))/kBT 
1 a^5s(/.i) 



a^SS](/.t) 

^j:{fi)/kBT 



3e \ii/kBT 



^Y.{^i)/kBT 



(16) 



Since 



cgY,{tJ.)/kBT general is nonzero, the thermopower 
in the clean limit is susceptible to impurity interaction 
and does not show universal behavior for a given chemical 
potential, different from what occurs in the case of finite 
magnetic field which we show in the next section. The 
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FIG. 2: DC conductivity a^xie^ /4:h) (left) and thermopower 
Sxx{kB/e.) (right) at fesT = 10~^eV vs chemical potential 
fJ.{eV) for y = 10 (blue lines) and V = W'^eV (dashed purple 
lines). The conductivity at ^ = for dashed purple line is 
around /2-Kh. The dotted brown line on the right figure is 
the first term in Eg. pSf) which is the dominant thermopower 
in large /i. 



DC conductivity and diagonal thermopower for general 
impurity strength are computed numerically and shown 
in Fig.(I2|). Similar to the results of self energy in Fig.(IT]) 
the DC conductivity shows asymmetry in fi for moderate 
impurity strength. The dip near /z = for V = IQeV 
case is the remnant signature of the clean sample, which 
shows up as a singularity in thermopower of the form 
Sxx oc for /i ~ 0. 



IV. FINITE MAGNETIC FIELD 

We assume the magnetic field B is perpendicular to 
the applied electric field E = Ex and choose the gauge 
A = (0, Bx sin 9, —Bx cos 9) . For E ^i) the eigenvalue of 
Hf) is given by the Landau level spectrum 



£;(7,n) 



-/i + "fvpl^^V^n + aB sm9 
-/i + 7e„ for n C N 



-/i 



gfisB sm9 = + €o (17) 



where a = g l^^/vpC^ 7 = ± is the band index, Ib — 
^Jl/eBsv[\9 is the magnetic length, and n denotes the 



index of Landau level. For Dirac electrons on the sur- 
face of HgTe a ~ 10~^/T. Thus for usual magnetic field 
strength we can in general ignore the Zeeman contribu- 
tion in the spectrum-'^l. Nevertheless we keep this Zeeman 
contribution in our calculation and the dominant feature 
that results from the Zeeman term is the lack of particle 
hole degeneracy at n = level. The loss of symmetry 
leads to anomalous features in off diagonal conductiv- 
ity (Jxy as well as thermopower. Within linear response 
we take i? ~ in the Hamiltonian and the field opera- 
tor ip{r) are expressed in linear combinations of energy 
eigenstates. 



V'(r) 



k 



i{k-\-g^B B cos B)y 



i{k-\-gfiB B cos 0)y / 







(x) 



with fy 



I exp 



g^isB sm9)/{"/vpy/2n/lB 

1 



-1 
2 



and Hn{x) is the Hermite polynomials^'- The sum over 
level index n is cut off by the linear spectrum boundary 
E{-i,n,nax) = -/i + iD. For D - 3.5eF, wf 5 x 10^ 
in the typical surface state of topological insulator— the 
upper cutoff obtained this way is rimax — 10^. In our 
numerical computation we take rimax = 2000 for prac- 
tical computation^® and checked that the results do not 
change much by comparing with rimax — 3000. In this 
representation is diagonal and the Green's function 
in energy eigenstate is 



Go(fc, n,7;iw) 



- E{-i,n) 

The impurity Hamiltonian, Eq.(l2]), in this basis is 



(19) 



Hi, 



V 



i—l p,k 



,i{k-p)yi 



i'*p,o{Xi)4>kfi{Xt)cl^oCkfi + ^ 



7/7 



r 



1 + fl 



7.7A/(l + A%.)(l + /7%n) 



(20) 



To include the effect of dilute random impurities we use 
the self consistent Born approximation and the averaging 
procedure over impurities positions as in Ref. m The 



position averaged Green's function is 
(G(fc,n,7;icj;ri,f2,f3, . . . ,rjv,„J) 
dr, 



[n 



1=1 "' 
G{k, n, 7; iuj) 



G{k, n, 7; iw; fi, ra, ^3, . . . , rjv™^) 



5 



In the averaging procedure we first integrate yi and then 
Xi and use the following identity 



dx(j)l ^{x)(f)k,m{x) = 5„ 



and then sum over momentum index in the propagator 
by using 



The self energy within self consistent Born approximation 



(21) 



In Eg. ([21]) the self energy, identity matrix 1, and Green's 
hmctions are all of dimension (2n„ + 1) x {2nm + 1) with 
n„i = rimax as the maximal Landau level cutoff by the 
linear spectrum. Denote 



fj^nf—'y.n) 



7,— 7,n 



(i + A',J(i + /^,„) 



and i„xn as n by n matrix with matrix element 1 along 
{rij, — JT-j}, the bare self energy term written in the energy 
eigenvalue basis is then given by 



Vlr 



Eo 









V 




-\-,n ^rij, 



Vlr 



An important point to note is that off diagonal term ap- 
pears when the Zeeman effect is important. Such terms 
are absent in graphene and are responsible for g-factor 
dependent contributions to transport. By performing an- 
alytical continuation {iuj — J- + iO^), the diagonal and 
off diagonal Green's functions (the off diagonal part de- 
scribes the change from 7 to —7 state or particle hole 
transition) are 



G(fc,n,7;w + 0+) 



w + ^ + 7e„ - E(7,n,a;) 



{uj + ^1- E(7, n, w) - E(7, -7, n, w))(w + ^ - E(7, n, uj) + E(7, -7, n, w)) - el 
Gik,0;oj + 0+) = — ^ , (22) 



UJ + ji — eo — E(0, (jj 
G(fc,n,7,-7;cj + 0+) = 



E(7,-7,n,w) 



(w + /i - E(7, n, uj) - E(7, -7, n, uj)){uj + fi - E(7, n, uj) + E(7, -7, n, uj)) - el 

I 



with diagonal self energy E(7, n, lo) at Landau level n, self 
energy E(0,cj) at zeroth Landau level, and off diagonal 
self energy E(7, — 7, n, oj) given by 



E(7, n, uj) = Vni/{1 - VZ{'j, n, io)) 
^{0,Lj)^VnJ{l-VZ{0,Lj)) 



(23) 



E(7, -7,n,w) = 



1 - Vj-j^nZ{-/, -7, n, uj) 



The explicit expression for internal propagators, 
Z{'y,n,uj), Z{0,uj), and ^(7,— 7,n,w) are given in the 
Appendix [C] Within these internal propagator the sum 
over k (in y direction) gives a factor Qc = V/2Trlg = 



2/{k^l%) = 2vj,/{D^l%), which accounts for de generacy 
of a given Landau level per unit cell and we drop the 
k dependence in the Green function in the Appendix [Cl 
The self consistent solutions of Eq. ip^ and Ea.([^ gives 
the density of states and the renormalization of Landau 
levels due to impurities^. In the case of weak magnetic 
field where the Zeeman contribution can be neglected we 
obtain the same results, modulo a constant factor due to 
spin and valley degeneracies in graphene, as in Ref. [itI . 



Our main goal is to study the transport properties of 
surface states. The particle current current correlation 
Lll.{uj) and Ll.y{uj) from Eq.© are given in the Ap- 
pendix|D] Before we proceed to evaluate the conductivity 
aij and thermal power Sij numerically for general impu- 
rity strength, we take the clean system limit F — ?> and 
study analytically the thermopower in the two dimen- 
sional helical metal. We show that in this clean limit we 
obtain universal feature for diagonal thermopower, simi- 
lar to its counterpart in the conventional two dimensional 
electron gas. 



A. Analytic results for weak impurities 

For weak impurities 1^ ~ we have 5G(fc,n;e) ~ 
7r(5(e + /x — e„). Taking this expression into Eq. (jDl[) gives 



rii I — 7"ii| 



= 



since the energy difference between Landau levels is large 
enough to prevent the transitions between different lev- 
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els without the help of impurity broadening. Also from 
with Ll\ given by 

de , , e , r 77r(5(e + /i — eo) 



4nllojh 



^tanh(^)^[- 



- /i + 70; — aei 



77r(5(e + /i — aei) 
e + ^ + 70; — eo 



E 



7:d{e + fj,- Ae„)7/2 
e + ^ + 7W - ae„+i 



TTS{e + 11- Ae„+i)7/2 



e + ^ + 7CJ - ae„ 



I tanh(- 



2T 



eo +0; 



eo - w 



(eo-a;)2-efJ (ei - eo) 



tanh( 



2T ' 



tanh(- 



2T ^ 



(ei + eo)^ - 



+ 



E- 
^ 2 



- (tanh( 



2T ' 



en + w 



tanh(^^)) 
e„ - uj 



(tanh( 



£n+i — f^ , 
2T ' 
e„+i + UJ 



tanh( 



- wj^ - e 



n+1 



2T 

-l-UJ 



■)) 



1 



ri+l 



)} 



(24) 



From Eg. dMl) it is easy to see that for eo = /i = we 
as a result of particle hole symmetry. For 



have Lll 
illL^o we have 



xy \'^- 



tj2tanh(^) 



— <^ - 2tanh( 



^ a;2tanh(^^) 



+(tanh( 



2T _ 

(ei+eo)2 



2T 

e„ 



2T 



tanh( 



2T 

2r 



(6l-eo)2 

^) _ tanh(^^^)) 



2T 



))]} 



(25) 



Here ujc = V2vf/Ib is the cyclotron frequency. Eq.((25|) 
suggests that the hall conductivity can be used as a probe 
for nonzero eo, or in other words, whether the Dirac spec- 
trum is gapped by the Zeeman term under the assump- 
tion that the sample is clean. By using Eg. ips]) in Eq. p^ 
we get 
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FIG. 3: Negative thermopower —Sxx as a function of chemi- 
cal potential /x for temperature ksT = 10~^ (Blue, Dashed), 
10"^ (Purple, Thick line), and 10"^ (Brown, Dot Dashed) in 
the clean helical metal. ^ and ksT are in unit of hwc and 
Sxx is in unit of kB/e. Around /i = —Sxx becomes or- 
der of 10^^ for all three cases. The peaks value are in nice 




agreement with Eq. (|26 
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FIG. 4: Negative thermopower ~Sxx as a function of chemical 
potential ^ for temperature fcsT — 10~^&Jc with gyromag- 
netic ratio g — 2 (blue thin line), g = 20 and B — lOT (purple 
thick line), and g — 20 and B — 20T (brown dot dashed line) 
in the clean helical metal, jj, and fcsT are in unit of fiuc and 
Sxx is in unit of ks/e. The peak height is not universal with 
large Zeeman effect. 



{[e(eo - e) - e(e - eo)](~2) - ( "^j^ + l)[e(ei - e) 



low temperature limit (i.e. huc 3> ksT) with chemical 
potential close to Landau level e„, is given by 



-e(e-ei)] + ( 



(ei + eo)' 



lei - eoj 
+ l)[e(e + ei)-e(-ei-e)] 



Tin 



^,,(/x,T^O) 



41n2^ 



if 



M = eo 



=p21n2^( 



-2 ^ [(e(e„ - e) - e(e - e„)) - (e(e„ + e) 

n=2 

-e(-e-e„))]} 
In the clean limit the off diagonal thermopower S^y 



3+( 



1^ 



if /i = ±ei 



^41n2ie. 



4W-2+( 

if = ie^v and iV > 2 



Syx = 0. The diagonal component Sxx 



Syy, in the 



(26) 

We can approximate (ei=Feo)^ — el{l + 2^) because 
q;(~ 10~^(l/r)) is small. Under this approximation we 
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see that for ^ = eq = g^BBsinO the thermopower be- 
comes — — In 2 — , which again serves as a measurable 
quantity for the Zeeman term. In real material, how- 
ever, this singular behavior around eo will be masked by 
the broadening of zeroth Landau level due to impurities 
or Coulomb interaction. For materials with negligible 
Zeeman interaction the thermopower reaches ^ for 

H = ie^v, irrespective of the strength of magnetic field. 
This is consistent with the Berry phase argumen t 
For weak magnetic field strength (eo <^ ei) the ther- 
mopower of two dimensional helical metal is similar to 
that of the conventional 2D metal in that there exists 
universality for peak value of thermopower in the clean 
limit^i^. The feature observed in the helical metal, how- 
ever, is different from the case of conventional 2D metal 
where every Landau level N has universal value — , . 

From Eq. (p6)) we see that for large Zeeman effect 
(ei — eo > iOc) the peak height of Sxx is no longer univer- 
sal and depends on the strength of Zeeman interaction. 
We artificially increases this Zeeman effect by increas- 
ing the gyromagnetic ratio and magnetic field strength. 
The results are shown in FiglH In addition to the non- 
universality another feature to note is the shift in the po- 
sition of the peaks refiecting the Zeeman contribution to 
the energy. Such a shift is most prominent for the lowest 
Landau levels, increasing the spacing for positive N. The 
significant reduction in the peak values as a function of 
magnetic field, as compared to graphene, provides a clear 
signature of the surface state. One consequence is that 
for the same impurity broadening, systems with larger 
Zeeman coupling possess gaps in the spectrum near the 
zeroeth Landau level even with the broadening from im- 
purity interaction. The effect of impurities are explored 
in the next subsection. 



B. General impurity strength 

For general impurity strength and large magnetic field 
strength we expand /-,^„ ~ (1 + ^^^f^ + 

Within this approximation -^,„ = "'^"^"'"^ The 
off-diagonal self energy S(7, —7, n, uj) is 

^ -VaB sine 

S](7,-7,n,a;)^ ^ ,.y„;.in/" i (27) 

The -irimax + 1 self energy terms in self consistency equa- 
tions Eq. ([22l) and Eq.([23|) are therefore simplified to two 
variables l](0,a;) = SoC"^) smd S(7, n,a;) = which 
we evaluated numerically in this section. We put band 
energy D = 3.5eV, Fermi velocity vp = 5 x lO^m/s, and 
impurity density n^mp — 10~^ in the following numerical 
computations. 

The real and imaginary part of the self energy are 
shown in FiglS] for B — 14T. For strong impurity 
strength (blue lines) the imaginary part of the self en- 
ergy shows particle hole symmetry (^ = in all self en- 



Re[E(6)]/fiu, Rem£)|/Jiu, 




FIG. 5: Self energy Eo(e) (left) and E(e) (right) vs e. Energy 
is in unit of fkoc and B — 14T in all figures. Top panel shows 
real part and bottom panel shows imaginary part. Impurity 
strength V = lO^hujc (blue lines), V — lO^&Jc (purple dashed 
lines), and impurity concentration Ui = 10~^. 




FIG. 6: Density of state per unit cell in unit of 1/hwc vs 
ui/ujc for B = 14r (left) and B = 7T (right). Blue lines are 
V — lO^hujc and purple dashed lines are for V = lO^fiuc. The 
brown dotted lines indicate the density of state with degen- 
eracy factor (jc and without impurities. 

ergy figures) while for moderate impurity strength (pur- 
ple dashed lines) the self energy shows more weight on 
the hole side. In clean limit, as discussed in last section, 
or with weak impurity strength the particle hole symme- 
try of the self energy is restored. This can be understood 
by taking — )■ 00 for vacancies and y for weak 
impurity in Eg. ipTj) . For F — >• cx), EoG > 1 and thus 
S = — n^G"^, which renders the self energy to be parti- 
cle hole symmetric. For y 0, E = niEo(l — SqG)^^ ~ 
n^So + j^jSqG. The first constant term can be absorbed 
in the shift of chemical potential and therefore the self en- 
ergy obtained by second order self consistency equation 
is particle hole symmetric. For intermediate impurity 
strength, the factor 1 in numerator of Eq. (|2ip is compa- 
rable with EqG. Thus for repulsive V the self energy is 
shifted to the left. 

From the self energy we obtain the density of state per 
unit cell as 

p(w) = Z5£cj[G(0; + 0+) + y G(n, 7; + 0+)] (28) 

TT — ' 

7,n 

The numerical results are shown in FiglHlfor two different 
magnetic field strengths. For V = IQphujc we see clear 
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FIG. 7: DC conductivity a^x (left) and a^y (right) as a func- 
tion of chemical potential ^. Conductivity is in unit of /h-n. 
H, ksT — 10~®, impurity strength V = 10^ (blue lines), and 
V = 10^ (purple dashed lines) are in unit of hujc- Magnetic 
field strength B = 14T for upper and B = 7T for lower fig- 
ures. The inset of upper right figure shows the blowup plot 
with smaller chemical potential range. The a^y in B = 14T, 
case has been shifted down by 50e^//i and a^y in 
B = 7T,V = lO^'&Jc case has been shifted down by 3.3e^//i 
to match the correct charge neutral point. 



signature of particle hole asymmetry and there are more 
holes than electrons when /i = 0. The Hall conductivity 
axy obtained in this case does not flip its signs at /x = 0. 
Rather the charge neutrality point has been shifted to 
positive II. Another feature is that the oscillation due to 
applied magnetic field is more apparent in > region 
than in a; < region for V = lO'^huJc- This feature also 
shows up in the longitudinal conductivity axx v.s. chemi- 
cal potential. The numerical results for DC conductivity 
ax'x and axy as a function of chemical potential fj, are 
plotted in FiglTl As is expected from the results of self 
energy and density of state, both axx and axy show more 
oscillatory behavior/larger steps for fi > region. 

The thermopower Sxx and Sxy obtained from Eq. (jl3p 
are shown in FigH] for B = 14T and B = 7T with 
V = lO^hujc and V = lO'^huc- The general features 
for Sxx are the peak positions corresponding to different 
Landau levels are shifted away from the n = Landau 
level and the width of the peaks increases due to the 
impurity broadening, similar to effect of temperature. 
There is no universal value at the peaks of Sxx due to 
impurities as in the case of conventional two dimensional 
electron gas^^. Near the charge neutrality point Sxy 
shows a peak while Sxx shows singular behavior around 
this point. The oscillatory behavior in Sxx and Sxy is 
more apparent for larger magnetic field, while the asym- 
metry is more prominent at lower fields. Note that the 
oscillations are entirely absent for negative chemical po- 
tentials for moderate impurity potentials. These results 
are in qualitative agreement with available scanning tun- 
neling microscopy(STM) data which find oscillations in 
density of states only for positive gate voltages^i^. 



Currently only limited experimental results for the 
thermopower of helical metals are available. Since our 
formulation is very similar to the case of graphene, we 
compare our thermopower results with data on graphene. 
In Ref. ,15 the observed peaks of magneto-thermopower 
Sxx shows the quantized l/N {N is the nonzero Landau 
level index) trend but with a reduced factor in height and 
peak broadening, which are qualitatively consistent with 
our numerically results shown in Fig|8]for V = lO^fiwc. 
In Ref. [l^ similar results have been reported, with the 
measured peak value of Sxx a.t N — —1 and B = 9T 
is AlfiV/K, which is smaller than the clean limit where 
Sxx = \n2kB/e = 0.69fcs/e ~ 59.7 ^iV/K. The ther- 
mopower shown in FiglS] for B = 14T, V — lO^hujc is 
similar to the case seen in Ref. [l^ that Sxx — 'iOiJ.V/K 
for = — 1 and smaller Sxx at |A^| > 1, suggesting the 
influence of impurity interactions. In both experiments 
a large peak at A^ = Landau level is observed which is 
consistent with our numerical results. This peak in Sxy 
can be explained by the presence of impurities in addition 
to the Dirac spectrum. 

The major difference between our results and experi- 
mental results in graphen o^^i^^'^'^ is the behavior of Sxx 
near A^ = region. Our numerical results suggest a large 
electron/hole peak between N = and A^ = 1/ — 1, show- 
ing remnant singular behavior around N = landau level 
in the clean limit as seen in Fig|31 In experimental data a 
hole/electron peak between A^ = and A^ = 1/ — 1 is ob- 
served instead. This discrepancy is related to the higher 
axx value when fi is at N — Landau level, as confirmed 
by the following observations. In all our numerical results 
the diagonal conductivity axx is close to the universal 
value e^/27r/i at /x = 0. By adding a Lorentzian profile 
to artificially increase axx around /x = we can repro- 
duce this results around A^ = seen in the experiments. 
Similar observations were reported in Ref. [2^ where they 
performed numerical study on discrete lattice model for 
graphene system. The physical reason for the increase of 
axx near A^ = Landau level is still an open question 
and we discuss the possible mechanism for this enhance- 
ment of axx near the Dirac node under magnetic field in 
Sec. El 



V. CONCLUSION 

In this article we compute linear response conductiv- 
ity and thermopower of the helical metal, both with and 
without magnetic field, in the presence of static, short 
ranged, random impurities. Analytic results of ther- 
mopower in the clean limit, i.e. impurity strength V ^ 0, 
are given in both cases. In the case of zero magnetic field 
the thermopower obtained is still susceptible to impu- 
rity interaction and does not show universal feature. In 
the case of finite magnetic field the thermopower shows 
universal value, independent of temperature and mag- 
netic field strength, when the chemical potential is fixed 
around Landau levels at low temperature. We also com- 
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fi/haic ft I tnii^ 



FIG. 8: Thermopower S^^ and S^y as a function of chemical potential /i. The thermal power is in unit of fc^/e and ^, fc^T, 
impurity strength V are in unit of hwc- Magnetic field strength B = 14 for upper two and B = 7T for lower two figures. The 
two left figures show the case of V = lO^huJc and two right figures show the case oi V = lO^fitJc. ksT = O.OS/kJc in all cases. 



pute the thermopower and conductivity numerically for 
the case of finite impurity strength. The effect of impuri- 
ties broaden the peaks in diagonal thermopower and the 
height of the peaks show no universal behavior, which is 
similar to the case in the conventional 2D electron gas^^. 
We also vary the impurity potential and study its affect 
on transport. For moderate impurity strength V we find 
the particle hole symmetry is broken while for large or 
small V the symmetry is restored in our perturbation 
formulation. 

Due to the lack of available data on thermopower mea- 
surement of the helical metal, we compare our ther- 
mopower results with experimental data on graphene. 
Our computation is on the helical metal system which dif- 
fers from graphene^ ^ mainly in the Zeeman effect. Since 
the Zeeman contribution is small for normal magnetic 
field strength (due to a ~ 10~'*/T) in both graphene 
and helical metal samples such as HgTe, our results are 
applicable to the case of graphene with short range, ran- 
domly positioned impurities modulo a numerical factor 
in conductivity due to valley and spin degeneracy in 
graphene. We find qualitative agreement between our nu- 
merical results of thermopower and experimental results 
in graphene, except for the Sxx near the zeroth Landau 
level. The discrepancy is due to the enhanced diagonal 
conductivity nearby zeroth Landau level. The enhance- 
ment in a^x seems to decrease with increasing magnetic 
ficld^^ . Possible mechanism includes the metal to insula- 
tor transition on the zeroth Landau level^l or electron 
hole puddles^*. Transport measurements of graphene 
near the minimum conductivity in Ref. [2^ partly sup- 
ports the latter mechanism but further study is needed 
to clarify this issue. Whether similar behavior is repro- 
duced in helical metals is an open question. 

The magneto-electric coupling, unique to topological 



insulators, leads to novel features in transport. The 
dependence of the thermopower on the Zeeman cou- 
pling provides a diagnostic of the state. To amplify the 
expected behavior we have analyzed thermopower and 
Nernst for materials with an artificially enhanced gyro- 
magnetic ratio. The relative importance of the Zeeman 
term to the Landau level splitting can be enhanced by an 
in-plane electric field, providing a possible route to de- 
tecting the signatures proposed. One caveat the electric 
field leads to additional effects such as non-linear trans- 
port, a detailed analysis of which is the subject of future 
investigation. 
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Appendix A: Derivation of self energy in zero 
magnetic field case 

The free electron Green's functions in the energy eigen- 
state basis Eq.® are 



SCBAll^, is 

(G's,s(fc,a;„;fi,r2,r3, . . . ,rAr.„„)) = G(fc,a;„) 



y 1/2 

^ iuj„ + /I - "fvpk 

iu!n + fJ- — JVpk 



= ]4iv— [ n J 



drj 



Gs^s{k,UJn;ri,'r2,r3, . . . ,rjVi„p) 



Thus the bare propagator spin-diagonal component) in 
Eq.Q is replaced by full propagator which leads to full 
self consistent Born approximationii 



7=±1 

Gi,l{k,UJn) = Gf^-f{k,UJn) 

In this basis, the impurity interaction in Eq.(l2]) is 



1 - Iv Efc,T,Nn + M - I'^pk - E(cj„)] 



-(Al) 



where S(aj„) is the electron self energy. This is inde- 
pendent of momenta due to spherical symmetry in the 
impurity potential in Eq.®. By using analytic continu- 
ation iLOn — > w -|- iO"'" we have the self energy of electron 
given by 



2A 



fe,7 



The impurity interaction is diagonal in the spin space. 
The position averaged Green's function, obtained within 



J 



Appendix B: Explicit form of n[e,V,nimp) 

The explicit form of K(e, V, riimp) after the momentum integral is given as 

K(e, y",n„„p) = (Bl) 
2ahD-^ (a2 - foS _ ^2) ^ (^2 ^ ^2^ _ ^)2 ^ 52) ^ ^)2 ^ ^2) ^^^^-1 ^a-D^ ^ ^^^^-1 ^o^D^ - 2 tan"! (f )) 



2ab + 2a2(fe - D)(& + Z?) + (62 + D^f^ 



with a = e - 5RE(e) and b = 5S(e). 



Appendix C: Expression for internal propagator terms 

The internal propagators for self energy computation, Z{'-^,n,uj), Z(0,a;), and Z{'y,—j,n^uj) are given by 
Z{j,n,uj) = 



r (1 + f'flm)9cG{m, 7; + 0+) ^ (1 - flJ^,mf-y,,n)9cGim, 7, -7; -I- 0+) ^ 



m,7 



(1 + + A,™) (1 + /2,J^(l + /2„.)(l + /2-_,J 



J 7,n 
1 + f2 



(7cG(0;tJ-H0+) 
/I 



_ 1. I _y -V m 



m,7 



gcfy,mf~y,mG{m, 7, -7; + 0+) J 
(1-f /2„)(l + /2. J 



+ 5cG(0;u; + 0+) 



(CI) 
(C2) 



s >r^r(l- /7,"/-7,"/^,m)5cG(m,7;a; + 0+) (1 + /^.„/_^.„/^,„/_^,,„)5cG(to, 7, -7; w + 0+) , 



%m (l + /|,™)^(l + /^V)(l + /-7,n) ^(1 + /2 „)(1 + /2^,„)(1 + /|,™)(1 + fl^^J 



f-y.n f—'y.n 



g,G(0;c^ + 0+) 



(C3) 
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Appendix D: Expression for L^^ and LJ.^ 



For the case of finite magnetic field in the energy eigenstates basis the particle current operators, as in Eq.®, are 



7/7," 



7/7 



'fc,n+l,7''fc,">7 



^ 

:.n,7>7 •\/(-'- + f'^,n){^ + /7,n+l) 



'^I,n,7'^fc,"+1.7 



} 



^7/7," 



(4,0Cfc,l,7 - 4,1,7'='=.0) + X! 



V(i + /7:«+i)(i + /7%) 



fc,n,7,7 y (1 + /7,n)(l + f^,n+l) 



} 



The particle current current correlation and Ll.\^{uj) from Eq.© are given by 



Ti? 1 /"^ He 

^ ' K(e + a;)-nf(£)]| y ^ [5^(0; e + 0+)SG(l, 7; e + + 0+) 



-y "^■'7,1 

^ 7,n 



+3G(1, 7; e + 0+)5G(0; e + + 0+)] + ^ ^^^^:2^i^-^3G(n + 1, 7; e + 0+) 

„ ^ , U + /7,n+in-'- + /7,nJ 



3G(n, 7; e + 0+)3G(n + 1, 7; e + w + 0+) 



E 

7 

E 



3G(0; e + 0+)3G(l, 7, -7; e + w + 0+) + 3G(1, 7, -7; e + 0+)3G(0; e + w + 0+) 



J 7,71 



n,7,7 V,-L -r- 



(l + /7%)A/(l + /7Vl)(l + /-7,n+l) 



V/(1 + /7,1)(1 + /-7.1) 

=5G(n + 1, 7, -7; e + 0+)9G(n, 7; e + w + 0+) 



E 



fj,nf—j,n 



^7 (1 + + + /-7,«) 



3G(n + 1, 7; e + 0+)5G(n, 7, -7; e + + 0+) 



E 



^7 (1 + /2„)^(l + /^Vl)(l + /-7.«+l) 
f^,nf—^,n 



3G(7i, 7; e + 0+)3G(n + 1, 7, -7; e + + 0+) 



E 



^7 (1 + /|,„+i)^(l + + 



3G(n, 7, -7; e + 0+)3G(n + 1, 7; e + + 0+) 



(Dl) 
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r |^tanh(-i;)[ ^ ^^[5JG(l,7;e + aa; + 0+pG(0;e + 0+) 



af2 

^ 7,n 



-9G(l,7;e + 0+)5iG(0;e + aa; + 0+)]+ ^ ^ '^'^^"^ ^ , jRG(n + 1, 7; e + ao; + 0+) 



G(n,7;e + 0+)- ^ ^ 



n,a=±l,7,7 ' •'l,n+l/\- • Jl,: 
•2 

SRG(n, 7; e + aw + 0+)S>G(n + 1, 7; e + 0+) 



a[3?G(l, 7, -7; e + aw + 0+)3G(0; e + 0+) - 5G(1, 7, -7; e + 0+)3fiG(0; e + aw + 0+)] 



a=±l,7 



y^(l + /2j(l + /2^J 



af2 



n,a=± 



^,7.7 (1 + /|,J^(l + /7%+l)(l+/V+l) 



^G{n + 1, 7, -7; e + aw + 0+)9G(n, 7; e + 0+) 



V 

,a=±l,7,7 (1 + J(l + /|,J(1 + /2^,J 



5RG(n + 1, 7; e + aw + 0+)3G(n, 7, -7; e + 0+) 



"/7,. 



- E 

„,«=±1,7,7 (1 + A'n) + fln+l)a + I-^,n+l) 

^fj,nf—j,n 



=3fJG(n, 7; e + aw + 0+)3G(n + 1, 7, -7; e + 0+) 



+ E 



,a=±l,7,7 (1 + /|,„+l)v/(l + /7,n)(l+/-7,n) 



5RG(n, 7, -7; e + aw + 0+)3G(n + 1, 7; e + 0+) (D2) 



The sum over k (in y direction) results in the normalization ^factor NV being replaced by 2Trl%. 
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